Symmetry induced Dynamics in 
four-dimensional Models deriving from 
the van der Pol Equation 



Ricardo Lopez-Ruiz^ 

^Department of Computer Science, Faculty of Sciences, 
^Institute de Biocomputacion y Fisica de Sistemas Complejos, 
Universidad de Zaragoza, 50009 - Zaragoza (Spain) 



Abstract 

Different models of self-excited oscillators which are four- dimensional extensions 
of the van der Pol system are reported. Their symmetries are analyzed. Three of 
them were introduced to model the release of vortices behind circular cylinders with 
a possible transition from a symmetric to an antisymmetric Benard-von Karman 
vortex street. The fourth reported self-excited oscillator is a new model which 
implements the breaking of the inversion symmetry. It presents the phenomenon of 
second harmonic generation in a natural way. The parallelism with second harmonic 
generation in nonlinear optics is discussed. There is also a small region in the 
parameter space where the dynamics of this system is quasiperiodic or chaotic. 
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1 Introduction 



Self-sustained oscillations are found very often in mechanical, electrical, biological, 
chemical and ecological systems pQ. Nonlinearities are required in order to attain 
this kind of behavior. The system reaches an stable oscillatory dynamics as a conse- 
quence of an internal balance between amplification and dissipation. For instance, 
the oscillatory patterns in some chemical reactions [2], the beating of the heart [3], 
the oscillating coexistence between different living being species jl] , the Benard-von 
Karman vortex street in the wake of a cylinder [S], etc. are examples of continuous 
and discrete systems modeled as self-excited oscillators. 

In this work, five different models of self-excited oscillators with certain simi- 
larities are collected. First, the fundamental brick of all of them, the van der Pol 
equation, is recalled. Second, the two-dimensional van der Pol equation is extended 
to a four-dimensional model in the complex plane. Two other models deriving from 
the former complex van der Pol system are revisited. They were investigated as 
models for the shedding of vortex structures in the flow behind a cylinder. In a next 
section, a new model of two coupled self-excited oscillators implementing the cor- 
rect symmetries to get second-harmonic generation is introduced and its qualitative 
dynamical behavior is sketched. Finally, the conclusions are presented. 

2 Four-dimensional extensions of the van der Pol 
Model: symmetries and dynamics 

The van der Pol equation, which displays a wide range of behavior from weakly 
nonlinear to strongly nonlinear relaxation oscillations, is one of the simplest models 
presenting self-oscillations. It reads: 




(1) 
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where x(t) is recording the time evolution of a harmonic oscillator perturbed by the 
nonlinear dissipative term e(l — x 2 )x. 

The system (P) presents the following symmetries: 

(i) Temporal translation symmetry: if x(t) is a solution of equation then 
x{t + 0), with a real constant, is also a solution of system ; 

(ii) Spatial inversion symmetry: if x(t) is a solution of the van der Pol equation 
then — x(t) verifies also the equation ((Tj). 

(iii) e-Parameter inversion symmetry: if x(t) is a solution of equation (PJ for the 
constant e then —x(—t) is also a solution for the parameter — e. Therefore, it is 
sufficient to find the solutions of equation JIJ for e > 0. The dynamics for e < is 
derived by applying this symmetry to the latter solutions. 

The dynamics of the van der Pol system is quite simple. It presents a fixed 
point at the origin and a unique periodic orbit for the whole range of the parameter 
e. For e < 0, this limit cycle is unstable, and, for e > 0, the limit cycle becomes 
attractive. Due to the uniqueness of this periodic solution, x(t), all the solutions 
generated from x(t) by the symmetries (i) and (ii) fall on the same orbit in the phase 
plane representation given by the coordinates (x(t),x(t)). Therefore, there exists 
a satisfying the equation x(t + 0) = — x(t), or equivalently, x(t + 20) = x(t). If 
we call T the period of the limit cycle, that is x(t + T) = x(t), the symmetries of 
the equation (P) impose the condition 20 = T. Then, the final dynamical regime 
verifies xit + T/2) = — x(t). Let us remark that this property is a consequence of 
the symmetries and the uniqueness of the stable periodic solution of system (0). 

The natural extension of the van der Pol equation to a four- dimensional model [Hj 
can be achieved by introducing in the equation (0) the complex variable z = Xi+ix2- 
Thus we obtain the complex differential equation 

z-e{l - \z\ 2 )z + z = 0, (2) 

where the real variables (xi,X2) are the coordinates of two symmetric coupled van 
der Pol self-excited oscillators which verify: ii — e(l — [x\ + x?,))xi + Xi = 0, with 
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i = 1,2. 



Equation (j2J) has the symmetries: 

(i) Temporal translation symmetry: z(t) — > z(t + (p), with a real constant; 

(ii) Phase symmetry: z(t) — > e l ^z(t) with ?/> real constant. It includes the spatial 
inversion symmetry: z(t) — * — z(£) when if) = 7r; 

(iii) Conjugation symmetry: z(t) — > with = — ^(t)). 

(iv) e-Parameter inversion symmetry: (z(t),e) — > (— z(— £), — e). 

The dynamics of system (j2J is damped to zero when e < and gives self-sustained 
oscillations when e > 0. These stable periodic motions present the form z(t) = 
e^r(t) with <p being a real constant and r(t) satisfying r — e(l — r 2 )r + r = 0. 
It can be represented as a straight line through the origin in the plane (xi,^) a ^ 
constant angle 0, r(t) oscillating along this line. Evidently, this solution verifies 
zit + T/2) = —z(t), with T the period of the oscillation. It is found that the trivial 
periodic solutions z = e ±lt are unstable. 

Equation (j2J) has been used in Ref. [H] as a fundamental brick to build a simple 
model for the release of vortices behind circular cylinders. The parameter e would 
depend on the Reynolds number of the flow: e < for Re < Re c and e > for 
Re > Re c where Re c represents the critical Reynolds number for the onset of the 
vortex shedding [HI El IB]- By introducing a phase symmetry breaking term controlled 
by the small parameter j3, 

z-e{\- \z\ 2 )z + z + i(3z = 0, (3) 

the dynamical properties of equation (J2J) are modified. Now there are only two 
oscillating solutions that asymptotically attract the flow: the symmetric mode 0i = 
(xi = X2) given by z\ = e n ^ri(t) and the antisymmetric mode O2 = (xi = —X2) 
given by Z2 = e -7r/,4 r 2 (t), where ri(t) and r 2 (t) verify r\— e(l— rf)ri+\\ — (— l) l /3]r.; = 
(z = 1,2). Both modes have also the property Zi(t + T/2) = —Zi(t), % = 1,2, with 
T being the period of the motion. If (5 > 0, Gi is stable and 02 unstable, then 
the oscillators (xi,X2) are in phase and the representation of a symmetric vortex 
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street is obtained j^j. If /3 < 0, Oi is unstable and 9 2 stable, then the oscillators 
(xi, X2) are out of phase and the representation of an antisymmetric vortex street is 
therefore obtained [7J|B|. 

A slight modification of the nonlinear coupling term in the dissipative force of 
equation Q produces the model: 

xi — e(l — x\ — (1 + j)xl))xx + x\ + fas = 0, 

x 2 - e(l - x\ - (1 + i)x\))x\ + x 2 + f3xx = 0, (4) 

where the parameter 7 controls the breaking of the conjugation symmetry. The 
residual symmetries of this system in the plane (xi,x 2 ) are: 

(i) Temporal translation symmetry: (xx(t), x 2 (t)) —>■ (x\(t + 0), x 2 (t + 0)), with 
<p a real constant; 

(ii) Diagonal reflection symmetry: (xi(t), x 2 {t)) — > (x2(t),xi(t)); 

(iii) (3-Parameter inversion symmetry: (xi(t), x 2 (t), (3) — > (x 2 (t), — Xi(t), — (3). 

A more general and 'robust' unfolding for the transition between the modes 0i 
and 02 is obtained in this new model [5 . Two different scenarios (Scenario I and 
Scenario II) for this transition are found. The nonlinear transition is performed by 
the mediation of two intermediary mixed modes: in Scenario I these mixed modes 
are unstable and in Scenario II they are stable. The Scenario II has been probably 
confirmed in the experiment performed by Zhang et al. with two flexible filaments 
interacting and oscillating in a flowing soap [Oj. 

3 The Model without Inversion Symmetry 

The inversion symmetry is a fundamental property of the solutions of the former 
models. We introduce here an inversion symmetry breaking term in the equation 
(j2J), namely the quadratic term (3z 2 . We sketch the drastic consequences of this term 
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on the dynamical behavior of the system. The new equation is 



z-e{l- \z\ 2 )z + z + fiz 2 = 0, 



(5) 



where fi is the coupling constant. The symmetries preserved in this system are: 

(i) Temporal translation symmetry: z(t) — > z(t + </>), with a real constant; 

(ii) Conjugation symmetry: z(t) — ► z(t); 

(iii) fi-Parameter inversion symmetry: (z(t),fi) — > (—z(t),—fi). Then it is suf- 
ficient to study the solutions of equation (J3J) for fi > 0. The dynamics for fi < is 
derived by applying this symmetry to the solutions found for fi positive. 

If we expand the equation (jHJ) in the real variables (xx, x 2 ), we have: 



When fi is different from zero, a perturbative analysis similar to that performed in 
Ref. [Hj for the model (J3J) does not give in this case any relevant information. We 
have performed the numerical integration of these equations and, roughly, two dif- 
ferent dynamical regimes are found. The regions where these dynamical behaviors 
take place are bounded by two curves, /?i(e) and fi 2 {^), in the parameter space (e, 0), 
with f}\ < fii- For instance, when e = 1, fix — 0.9 and fii ~ 1.97, and the general 
characteristics of these regions are: 

(a) Region I : if < fi < fix or fi > fi 2 , a unique stable periodic orbit exists in 
phase space. It is found that if the coordinate x 2 oscillates with period T 2 = T then 
the coordinate x\ oscillates with the half of a period T\ = T/2. The system presents 
in a natural way the phenomenon of second-harmonic generation in this parameter 
region (Fig. 1(a)); 

(b) Region II : if fix < fi < fi 2 , the dynamics of the system displays three dif- 
ferent behaviors in the four-dimensional phase space: quasiperiodicity (Fig. 1(b)), 
periodic windows where the relation T 2 = 27\ is also verified (Fig. 1(c)), and chaos 



Xx — e(l — {x\ + x\))xx + Xx + fi{x1 — x\) 
x 2 - e(l - (x\ + x\))x 2 + x 2 + 2fixxx 2 



0. 



(6) 



(7) 



(Fig. 1(d)). 
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The second-harmonic generation phenomenon in Region I is a consequence of two 
facts: the symmetries of the system and the existence of only one stable periodic 
orbit. Let T be the period of this orbit, then (xi(t + T), x 2 (t + T)) = (xi(t), x 2 (t)). 
Due to the uniqueness of the stable limit cycle, the solutions generated from the 
time translation and conjugation transformation of this periodic orbit draw the 
same path on the plane (xi,x 2 ). There exists therefore a verifying that (xi(t + 
4>) , x 2 (t + <ft)) = (x\(t), — x 2 (t)). As in the van der Pol system, the only coherent 
solution for the equation x 2 (t + (fi) = —x 2 {t) is (j) = T/2. Then, x 2 (t + T/2) = —x 2 (t) 
and x\(t + T/2) = x\(t), that is T 2 = 2T\ = T. The oscillator x\ performs two 
oscillations for each oscillation of the oscillator x 2 . The limit cycle on the plane 
(xi, x 2 ) presents hence a pair number of lobes in the direction of the x 2 axis and it 
shows an eight-folded shape (Fig. l(a-c)). 

Let z(t,/3) be the unique stable periodic orbit in the Region I. By numerical 
inspection we find out that z(t,/3 = 0) = z (t) = (0,X2,o(^)), with X2,o(£) verifying 
x 2,o — e (l — x \ 0)^2,0 + £2,0 = 0. In order to perform an alternative perturbation 
analysis [10] in the parameter /3, we expand z(t,/3) = (xi(t, (3), x 2 (t, (3)) in power 
series of (3: 

00 

Xl (t,f3) = E^x^), (8) 
i=i 

00 

x 2 (t,P) = x 2>0 {t) + y £l3lx 2>j (t), (9) 

j=i 

where Xij{t) (i — 1,2 ; j — 1, . . . , 00) is calculated by an iterative process as a 
function of x 2t0 (t). Then, x it j(t) = Xij(x 2t0 (t)) , with i = 1, 2 and j — 1, ... ,00. The 
tilde of Xij indicates the different functional dependence of Xij on x 2t o(t). Thus, we 
have: 

00 

xi(x 2 , (t),/3) = Y,P j xiA x 2,o(-t)), ( 10 ) 

00 

x 2 {x 2fl {t),(3) = x 2>0 {t) + Y i l3 i x 2>j (x 2fi (t)). (11) 

i=l 

The /3-parameter inversion symmetry (z(t),(3) — > (—z(t),—f3) implies that the 
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relations: 



xi(-x 2 , (t),-(3) = -x 1 (x 2 ,o(t),/3), (12) 
x 2 {-x 2>0 (t),-P) = -x 2 (x 2 ,o(t),P), (13) 

must be verified by the correct power series expansion. Furthermore, the conjugation 
symmetry (z(t),/3) — > (z(t),j3) imposes the conditions: 

xi(-X2,o(*),/9) = xi(x 2 ,o(t),/5), (14) 
x 2 (-x 2t0 (t),/3) = -x 2 (x 2)0 (t),(3). (15) 

Retaining the correct terms which verify the parity conditions (jl2H15|) . we arrive at 
the correct power series expansion: 

oo 

= E^'" 1 ^--^ (16) 

i=i 

oo 

#2 = ^2,0 + J2P 2j X 2 ,2j, (17) 
J'=l 

where £2,2^ are odd functions of x 2 o and Xipj-i are even functions of x 2 o for all j. This 
last property implies that x\ is a function of x\ . Because x 2 ${t + T/2) = — i 2 ,o(t) 
is verified we find that x| (^ + T/2) = x 20 (t). Therefore the frequency associated 
to the oscillator Xi(t) is the doubled frequency of the oscillator x 2 (t), as it has been 
shown above by the argument of the symmetries. 

As an example, the first terms of the expansion (jl6H17j) . x\ ~ and x 2 ~ 

x 2 ,o + P 2 x 2t2 , satisfy the differential equations: 

x 2 ,o - e(l - xl fi )x 2fi + x 2 ,o = 0, (18) 
x hl - e(l - xl fi )x 1A + x 1A = xl >0 , (19) 
x 2 , 2 - e(l - ^2,0)^2,2 + [1 + e(z! )0 )]x2,2 = -2x2,0^1,1-^2,0^1,1- ( 20 ) 

We have verified the correctness of this expansion numerically, by recording x 2i o(t) 
from equation (fT%|) and by reinjecting this signal in equations ()19H2()j) . thus obtaining 
Xi i(t) and x 2j2 (t). The comparison of this approximation with the real dynamics of 
x\(t) and x 2 (t) derived from equation (jHJ) presents a good agreement for (3 slightly 
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different from zero. Also, the dependence on the parameter f3 of the limit cycle 
period has been calculated by computational methods for different values of e. The 
results are shown in Fig. 2. 

The region of the parameter space where the system is chaotic is also investi- 
gated by computational methods. We apply directly the method of calculating the 
sensitive dependence on initial conditions for each particular pair (e, (3) of parameter 
values. When the system shows exponential dependence on the initial conditions, 
that is, when two infinitesimal separated initial conditions become separated of the 
system size order, then the system is chaotic and the pair of parameters is collected. 
The location of the chaotic region in the parameter space (e, 0) is shown in Fig. 3. 

Second-harmonic generation is one the first phenomena studied in nonlinear op- 
tics [H]. Its discovery [T2] put in evidence the possibility to modify the optical 
properties of a material system under the action of an applied optical field suffi- 
ciently intense (laser). Thus, in that process, a fraction of the power in the incident 
radiation at frequency w is converted to radiation at the second-harmonic frequency 
2w (Fig. 4). One of the causes of this phenomenon is the non- vanishing second-order 
susceptibility x °f the crystal on which the laser beam is incident. This possibility 
is only achieved in crystals which do not display inversion symmetry. Therefore one 
would be tempted to establish some formal relation between the behavior of the 
model (JHJ) in the Region I and that optical system (Fig. 4). From a qualitative 
point of view, this model presents the breaking of the inversion symmetry and the 
frequency doubling phenomenon in a natural way. 



4 Conclusions 



In the past few years, important geometrical and topological studies have been de- 
veloped for understanding the periodic orbit organization of low dimensional chaotic 
attractors [inilllj. However, we still experience, for instance, some difficulty to get 
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a good visual understanding of the entangled dynamics of a particle following a 
chaotic trajectory in a strange attractor, even if the system is a three dimensional 
one. Probably, this a consequence of the lack of some future technical development 
that would improve and would approach better this visual representation more than 
a consequence of the lack of theoretical tools for its comprehension. 

Thus, the study of four- dimensional systems can be a challenge in the next future. 
In this work, different self-excited oscillators which are four-dimensional extensions 
of the van der Pol equation have been reported. Symmetries play a determinant role 
in the dynamics of all of them. In the model (J2J), the circular symmetry forbids the 
capability of this system to select a final dynamical state with a fixed symmetry. 
By introducing a phase-symmetry breaking term, i/3z, and depending of the sign 
of the control parameter /3, the dynamics in the model (JHJ) can settle down on two 
different oscillation states: an in-phase mode 0i and another out-of- phase mode 02- 
Each mode might represent a symmetric and an antisymmetric Benard-von Karman 
vortex street, respectively j^J. In the model an additional breaking circular 
symmetry term in the dissipative force, monitored by the parameter 7, provokes 
the existence of two possible scenarios for the nonlinear transition between the two 
oscillation modes, Oi and 62. 

Let us observe that the final model © implements the necessary breaking term 
of the inversion symmetry, f3z 2 , in order to obtain second-harmonic generation when 
two identical self-excited oscillators interact. When (3 = the self-excited oscilla- 
tors are uncoupled and they oscillate to the same frequency w. It can represent a 
laser beam before reaching a crystal with non- vanishing second-order susceptibility. 
When laser beam passes through the crystal a nonlinear interaction takes place with 
the medium. Then the two self-excited oscillators are coupled (/? ^ 0). One of them 
preserves the original frequency w and the other one oscillates to the doubled fre- 
quency 2w (Fig. 4). Other more complicated phenomena such as quasiperiodicity 
and chaos are also found in a small and very localized region of the parameter space 
of the model (0) (Region II ). In the case of its experimental existence, the fact 
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that the two signals become chaotic in this regime as a consequence of the nonlinear 
coupling caused by the crystal would convert this unseen behavior in a surprising 
phenomenon. 
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Figure Captions 



Fig 1. Projections on the plane (xi,X2) of different stable attractors of system 
(jo]), (a) Projections of the stable limit cycle for e = 1 and (3 = 0.15. The biggest 
orbit represents the projection on the plane (£2,^2)- The smallest one is the pro- 
jection on the plane (xi,X\). The other eight-shaped orbit is the projection on the 
plane (xi,^). (b) Projection of the quasiperiodic attractor for e = 1 and (3 = 1. 
(c) Projection of a stable orbit in a periodic window for e = 1 and (3 = 1.2. Also, in 
this case, T2 = 27\. (d) Projection of the chaotic attractor for e = 1 and (3 = 1.6. 

Fig 2. (a) Plot of the normalised period, T/2tt, of the limit cycle of the van der 
Pol oscillator as a function of the nonlinear friction strength given by the parameter 
e. (b) Plot of the normalised period, T/2n, of the limit cycle of system (jSJ) as a 
function of the coupling parameter (3 for e = 0.5 (triangles), e = 1 (circles), e = 1.5 
(diamonds) and e = 2 (squares). 

Fig 3. Values of the parameters (e, (3) where the system (JSJ) is chaotic. The sen- 
sitive dependence on initial conditions has been verified by computational methods. 

Fig 4. Schematic representation of the second-harmonic generation phenomenon 
in nonlinear optics. 



14 



